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The acceleration of protons and electrons in a reconnecting current sheet (RCS) is simulated with 
a particle- in-cell (PIC) 2D3V code for the proton-to-electron mass ratio of 100. The electro-magnetic 
configuration forming the RCS incorporates all three components of the magnetic field (including the 
guiding field) and a drifted electric field. PIC simulations reveal that there is a polarisation electric 
field that appears during acceleration owing to a separation of electrons from protons towards the 
midplane of the RCS. If the plasma density is low, the polarisation field is weak and the particle 
trajectories in the PIC simulations are similar to those in the test particle (TP) approach. For the 
higher plasma density the polarisation field is stronger and it affects the trajectories of protons by 
increasing their orbits during acceleration. This field also leads to a less asymmetrical abundances 
of ejected protons towards the midplane in comparison with the TP approach. For a given magnetic 
topology electrons in PIC simulations are ejected to the same semispace as protons, contrary to 
the TP results. This happens because the polarisation field extends far beyond the thickness of 
a current sheet. This field decelerates the electrons, which are initially ejected into the semispace 
opposite to the protons, returns them back to the RCS, and, eventually, leads to the electron ejection 
into the same semispace as protons. Energy distribution of the ejected electrons is rather wide and 
single-peak, contrary to the two-peak narrow-energy distribution obtained in the TP approach. In 
the case of a strong guiding field, the mean energy of the ejected electrons is found to be smaller 
than it is predicted analytically and by the TP simulations. The beam of accelerated electrons is 
also found to generate turbulent electric field in a form of Langmuir waves. 



1. INTRODUCTION 

The observations in hard X-rays and gamma-rays of 
solar flares or geomagnetic tail imply that the essential 
part of the released energy has to be converted into ac- 
celerated particles - electrons and ions. It is usually as- 
sumed that the acceleration occurs as a result of conver- 
sion of the free magnetic energy released in a magnetic 
reconnection. The following three mechanisms of parti- 
cle acceleration are normally considered: acceleration by 
an electric DC field [l|, stochastic acceleration of differ- 
ent kinds 3 and shock acceleration [3)]. However, it is 
not still clear which of the mechanisms of particle accel- 
eration really takes place in solar flares or geomagnetic 
tail [4, 5]. Each mechanism can be found contributing 
to the emitted radiation and in order to distinguish their 
effects more accurate simulations of particle trajectories 
and energy spectra of accelerated particles are required. 

In the present study we focus on particle acceleration 
by electric fleld inside a reconnecting current sheet. This 
field can be considered as a logical step in progressing 
with a conversion of the magnetic field energy into the 
energy of accelerated particles. The first analytical study 
of this acceleration mechanism was proposed by Speiser 
@ who found that the particles can be accelerated to 
rather high energies after they enter the configuration 
with reconnecting magnetic fleld lines forming an RCS 
and a perpendicular (drifted) electric field. This theory 
was extended by Litvinenko and Somov [1,], Litvinenko 
01 who considered an RCS with all three components of 
the magnetic field and carried out analytical estimations 
for the energy of accelerated particles. In particular, the 
authors found that the guiding magnetic field can signif- 
icantly increase the energy of accelerated particles. 



Efthymiopoulos et al. [8| have analytically studied the 
trajectories of particles inside a similar current sheet con- 
figuration by means of the dynamical systems methods. 
In particular, they found that in a certain magnetic con- 
figuration a particle can be trapped inside the RCS. Con- 
ditions of a particle trapping have been found in case of 
a weak guiding magnetic field, while for a strong guiding 
field particles always follow escaping orbits. The energy 
gain for escaping particles is determined as a function of 
initial position and velocity. 

Particle trajectories in the similar 3D magnetic config- 
uration with the guiding field were further studied nu- 
merically by usinga test particle approach by Zharkova 
and Gordovskyy (constant electric field) and Wood 
and Neukirch (electric field enhanced near the X- 
nuUpoint due to anomalous resistivity). Zharkova and 
Gordovskyy [§] showed that the trajectories of particles 
with the opposite charges (electrons or protons) can be 
either fully symmetric or strongly asymmetric towards 
the midplane of the RCS depending on the ratio between 
the magnetic field components. If the guiding field is 
strong enough, the accelerated electrons and ions are 
found ejected into the opposite directions with respect 
to the midplane causing charge separation, if the guiding 
field is weak, the particle are ejected in even proportions 
into each side from the midplane. The direction of ejec- 
tion is shown to be dependent on the directions of the 
transverse and guiding magnetic fields as well as the di- 
rection of the drifted electric field. It was also shown that 
the accelerated particles gain energies up to 100 keV for 
the electrons and up to 1 MeV for the protons for the 
electric field magnitude of 100 V/m 0, El- 

As a result, the energy spectra of accelerated particles 
also depend on a magnetic field topology, an electric field 
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strength and the dependence of transverse magnetic field 
variations on a distance from the X-nuUpoint. For exam- 
ple, the spectral indices of energy spectra vary from 2 for 
electrons and 1.5 for protons if the transverse magnetic 
field linearly increases with the distance, or they change 
for a particular magnetic model adopted from the MHD 
simulation of Somov and Oreshina fl^ to 1.8-2.2 for elec- 
trons and 1.3-1.8 for the protons. In the similar model 
Wood and Neukirch [l^l obtained for electrons a power- 
law spectra with index 1.5. The all spectral indices for 
both electrons and protons become much higher if the 
transverse magnetic field increases exponentially with a 
distance from the X-nuUpoint (Tsj . 

Therefore, particle trajectories gained during accelera- 
tion in an RCS is usually considered non-self-consistently, 
i.e. in modelled electric and magnetic fields. Such stud- 
ies can give some information about accelerated parti- 
cles but they do not include fields generated by the ac- 
celerated particles. In the previously mentioned papers 
the current sheet was assumed to be already formed and 
remained stationary while reconnecting. On the other 
hand, there are also studies where magnetic reconnec- 
tion is treated in a framework of kinetic (self-consistent) 
simulation [see, e.g. [13, [H, [H, [i3|- Most of this studies 
are focused on the reconnection itself: determining the 
rate of reconnection and the outflow velocity, studying 
the structure of the dissipation region. 

However, only in recent studies the PIC simulations 
were applied to consider particle acceleration during a 
magnetic reconnection in the vicinity of X-nuUpoint. 
Tsiklauri and Haruki [l3| used PIC simulation to study 
the whole evolution of reconnecting magnetic field which 
started from the stressed X-nuUpoint configuration and 
evolved to the unstressed one within a few thousands of 
Alfven times. While the main goal was to investigate the 
reconnection processes in the vicinity of X-nuUpoint in a 
compressed current sheet, the authors show that some 
fraction of the the released magnetic energy is trans- 
formed into the kinetic energy of accelerated particles. 
The energy spectrum of accelerated electrons is power- 
law with the index in a range from 4 to 5.5 and maximal 
energy from 100 keV to 2 MeV depending on the stress 
factor of the X-nuUpoint. Another model of particle ac- 
celeration was also studied by using the PIC simulation 
by Drake et al. [lB| in the reconnecting magnetic con- 
figuration called magnetic islands. The PIC simulation 
reveals that the multiple interactions of electrons with 
these magnetic islands allow them to reach relativistic 
energies. However, there is no evidence that such mech- 
anism can be effective for ions. 

In order to study particle acceleration in a current 
sheet in the current paper we will also use the PIC ap- 
proach, which takes into account electric and magnetic 
fields generated by accelerated particles. Similar to the 
other authors, we do not simulate the whole reconnec- 
tion site but a relatively small part of the RCS at some 
distance from the X-nuUpoint with the size large enough 
to contain a whole particle orbit before its ejection. Sim- 



ilar to Drake et al. [lB|, we study the 3D current sheet 
which is already formed by a reconnection process. The 
RCS configuration in our study is similar to one discussed 
by Litvinenko and Somov f]\ , Zharkova and Gordovskyy 
[9| . Since the time scale of particle acceleration is much 
smaller than the typical MHD time scales, we assume 
that the background magnetic field which forms the RCS 
does not evolve during the particle acceleration. 

Although it has to be noted that PIC simulations for 
space plasmas have a few essential problems [see for the 
full details introduction inil8i]. Probably, the main prob- 
lem for the explicit PIC simulations is a limited size of 
the simulation region restricted by a modern computer 
power. This problem is the result of a condition that the 
step of the spatial grid cannot exceed the Debye length 
Ad — {kT/Anne^y^'^ . For example, in the solar corona 
the Debye length is of the order of 10^'^ m. Thus, the 
simulation regions in the PIC simulations of the mag- 
netic reconnection carried out by Birn et al. p^ . Tsik- 
lauri and Haruki [l3| in the solar corona conditions would 
have sizes of the order of few metres at the most. There- 
fore, it is very doubtful that PIC simulations can be used 
to study large scale reconnection events often occurring 
in the Sun Il9|. 

Much more realistic approach is to use the PIC sim- 
ulation to study small scale objects, with the sizes of 
the order of a current sheet thickness. For example such 
the study was carried out by Drake et al. [lH for the 
PIC simulation of two magnetic islands within a current 
sheet. However, the size problem discussed above leads to 
the following restrictions. For the current sheet thickness 
equal to the ion inertial length. Si, the number of cells 
across the current sheet in a PIC simulation have to be 
Si/)^D = Cy^mi/{kT), which is 3 • lO'^ for the solar corona 
temperature or 3 • 10^ for the magnctosphere. In order 
to reduce this number, Drake et al. 16] used a reduced 
magnitude for the speed of light, c — 20Va ~ 6 • 10^ m/s, 
where Va is the Alfven velocity in the magnctosphere. 
Another way to reduce the number of cells was used in 
the PIC simulation carried out by Karlicky [g^] , who con- 
sidered the high temperature electron-positron plasma, 
for which the ratio St/Xn was as low as 10. 

We assume that the current sheet is formed by the 
plasma with the typical coronal density of 10^^ m~'^, 
which sustain the equilibrium with the background mag- 
netic field. In order to avoid the problem with the small 
Debye length, only a small fraction of the plasma par- 
ticles (with density of 10^" m^'^) is included in the PIC 
simulation. This makes the ratio Si/Xd to be of the or- 
der of 10. This approach allows us to have a reasonable 
number of cells and super-particles without using the re- 
duced speed of light or ions with very small mass. Also, 
this approach can be considered as a development of the 
TP method towards the full particle simulation. As it 
is shown further, even such a small number of simulated 
particles can demonstrate the new effects, which are not 
possible in the TP approach, for example, they can gen- 
erate a strong polarisation electric field. We also use the 
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FIG. 1: Simulation model. 



TP approach to determine how the locally induced fields 
affect the particle trajectories. 



2. THE MODEL 

The model accepted in this paper (Fig. [1]) is a further 
development of the test particle model used by Zharkova 
and Gordovskyy [§]. We consider the reconnecting mag- 
netic field as a background field, which is caused by the 
external processes of magnetic reconnection and intend to 
simulate with the PIC approach the particle trajectories 
and their density/energy spectra in this field. Contrary 
to the TP simulations, plasma particles in the PIC simu- 
lations are considered to generate their own electric and 
magnetic fields, which is now self-consistently taken into 
account. 

TP simulations have shown that the acceleration time 
is of the order of 10^^ s for the electrons and 10~^ s for 
the protons. Since this time is much shorter than the 
time of the reconnecting magnetic field variation [2l| . 
then we can assume that during simulation the back- 
ground magnetic field is stationary. Also, from the TP 
simulations we conclude that travel distances of acceler- 
ating particles along the RCS are of the order of 10 km at 
most (for the protons) [l^ . Thus we can assume that this 
length is much shorter than the length scale of the mag- 
netic field variation along the current sheet. In addition, 
as it is generally accepted, we suppose that the magnetic 
field variations across the current sheet has much shorter 
length scale than its variation along the current sheet, 

L^<^L^,Ly. (2.1) 

Thus, our simulation domain is a small part of the 
RCS (see Fig. [T]), large enough to contain the full trajec- 
tories of accelerated particles. The background magnetic 
field is stationary and vary inside this domain only in 
the X direction, which is perpendicular to the RCS. We 
take into account all three components of the background 
magnetic field. The main component depends on x 



as follows: 

B,(x) = -B,otanh(^-^^ . (2.2) 

The Bx component is assumed [like \yS^ to be constant 
inside the simulation domain: 

= —Bxo- (2.3) 

The guiding (out-of-plane) magnetic field By is maximal 
in the midplane and vanishes outside the RCS: 

By{x) = ByQsech (^-^^ . (2.4) 

Note, that if ByQ — the configuration corresponds to 
the Harris sheet equilibrium, and if ByQ = B^o the equi- 
librium is force- free. 

The inflow of plasma into the RCS combined with the 
condition of the frozen-in magnetic field leads to the in- 
duction of the drifted (out-of-plane) electric field Ey. In 
order to provide the inflow of plasma in our simulation 
domain we set up a background electric field, as those 
drifted in with velocity Vi„ by a magnetic diffusion pro- 
cess [U. 

Ey = Eyo ~ BzoVin, (2.5) 

where Vin is the infiow velocity. 

In our study we use the following values for the cur- 
rent sheet parameters: the main component of the mag- 
netic field Bzo = 10~^ T, the current sheet half-thickness 
Lx ^ 1 m, the drifted electric field Eyo — 250 V/m and 
the guiding, Byo, and transverse, B^o, components of 
the magnetic field are varied to study how they inffuence 
particle acceleration. 

3. TEST PARTICLE SIMULATIONS 

For better understanding of particle trajectories inside 
an RCS during the PIC simulations let us first investigate 
the results obtained in the TP approach [similar to those 
byi. 

A typical trajectory of a plasma particle (Fig. ^ in 
electro-magnetic configuration shown in Fig. [1] consists of 
the three parts: E x B drift, acceleration by the electric 
field and ejection. (i)Outside the current sheet, where 
the magnetic field is strong, particle motion is adiabatic 
and can be described as the superposition of a magnetic 
gyration and a drift in the orthogonal electric and mag- 
netic fields towards the midplane of the RCS. (ii)Inside 
the current sheet the particle moves along the y axis and 
is accelerated by the electric field Ey. The exact tra- 
jectory of the particle depends on the value of guiding 
field By Q. For a large By the particle remains mag- 
netised and reaches higher energy than for a small By, 
when the particle becomes unmagnetised inside the cur- 
rent sheet. (iii)When the particle gains enough energy 
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FIG. 2: Trajectories of two protons (solid) and two 
electrons (dashed) on the x-V^ phase plane during 
acceleration inside a current sheet (TP simulations). 
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the magnetic configuration cannot confine the particle 
motion any more and it is ejected from the RCS. Outside 
the current sheet the particle motion is again adiabatic. 
However, due to a high velocity of the particle along the 
magnetic field line and non-zero the x component of 
the particle velocity is larger than the E x B drift velocity 
and the particle moves away from the RCS. 

Fig.[2]shows the trajectories in the x-Vz phase plane of 
two protons and two electrons (one proton and one elec- 
tron inflow from the a: > semispace another pair from 
the X < semispace). The mass ratio is chosen to be 
nip/me = 100 for the comparison with the PIC simula- 
tions in further sections. The magnitudes of the magnetic 
field components By and B^ are chosen such that elec- 
trons remain magnetised in the vicinity of the midplane 
while protons arc unmagnetised during the acceleration 
phase. In this case the energy of the unmagnetised pro- 
tons upon ejection can be estimated [see 7] as 



coming from a; > and electron coming from x < in 
Fig.©. 

As it is seen in Fig. O "transit" and "bounced" parti- 
cles gain different energy, and this effect is stronger for 
electrons. Note, that the difference in motion of "tran- 
sit" and "bounced" particles vanishes when By ^ [see 
also Figs. 6,7 inU^. On the other hand, if B^ = parti- 
cles cannot escape from the RCS and their acceleration 
is only limited by the size of the current sheet in the y 
direction. 

Since for the magnitude of By and B^, which were used 
in Fig. [2] the electrons are magnetised, then as it was 
shown by Litvinenko 01 they follow magnetic field lines 
and gain the energy: 



yO 



(3.2) 



However, this relation is only valid for the "transit" mag- 
netised electron. The "bounced" electron, on the other 
hand, experiences a repelling force caused by Ey while 
travelling along the magnetic field line towards the mid- 
plane. When this force overcomes the attraction due to 
the E X B drift the "bounced" electron is getting ejected. 
Since it cannot reach the midplane the "bounced" elec- 
tron gains less energy than the "transit" one. 

The ejection energies of electrons are plotted in Fig.[5al 
and I3bl as a function of B^q and Byo, respectively. The 
lines correspond to Eq. (|3.ip (dashed) and Eq. (|3.2p 
(solid) obtained analytically by Litvinenko in the lim- 
its of weak and strong By. The energies of "transit" 
electrons are in a good agreement with the analytical es- 
timations, i.e. with Eq. ()3.ip for weak By and B^ when 
electrons are unmagnetised, and with Eq. ()3.2p for strong 
By and B^ when electrons are magnetised. On the other 
hand, energies of the "bounced" electrons for any By and 
Bx coincide with these of the unmagnetised ones given 
by Eq. (|3.ip . This means that By enhances the ejection 
energy of the "transit" electrons only. 
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BxQ 



(3.1) 



where TOq is the particle mass. As it was shown by 
Zharkova and Gordovskyy Q if ByQ is strong enough 
(> 1.5 • 10~^Bzo) all protons (regardless of side they en- 
tered from) are ejected to one semispace while all elec- 
trons are ejected to the opposite semispace with respect 
to the midplane. Indeed, in Fig. [2] both protons are 
ejected to the semispace with negative x and both elec- 
trons to the a; > semispace. In order to distinguish 
two types of trajectories, the particles that inflow from 
and are ejected into the same semispace will be referred 
to as "bounced" particles (this corresponds to the pro- 
ton coming from a: < and electron coming from a: > 
in Fig. [2). Particles that are ejected into the opposite 
semispace from the one they come from will be referred 
to as "transit" particles (this corresponds to the proton 



4. PARTICLE-IN-CELL SIMULATIONS 

4.1. Problem formulation 

Although the TP simulations can provide some valu- 
able results of particle motion inside an RCS, it does not 
take into account the electric and magnetic fields gen- 
erated by accelerated particles in the simulation region. 
In order to include these fields, we used 2D3V particle- 
in-ccU simulation code developed by Verboncoeur and 
Gladd [22]. PIC method, similar to TP, is based on the 
equation of motion for plasma particles 



dt 



TO,, 



E -f E -I- V, 



B + B 



(4.1) 



where besides of the background fields E and B given 
by Eqs. (I2.2p - (|2.5p . the local self-consistent fields E and 



5 



10° 



10 



o 

4 

>■ 10 


10^ 



10 



4^ 
















X 














- 



































10 



10 



B,Tesla^° 



10 ■ 



(a) Byo = 10-4 T 




B , Tesia 

V 

(b) B^o = 2- 10-4 T 

FIG. 3: Electron energy in the TP simulations upon 
ejection as a function of {a.)Bxo and {h)Byo. Crosses - 
"transit" electron, circles - "bounced" electron, dashed 

line - Eq. ^J^, solid line - Eq. ([321). B^o ^ IQ-'"^ T, 
Sj^o = 250 V/m, = 1 m. 



B induced by the accelerated particles are taken into 
account. These fields are calculated from the Maxwell 
equations: 



— = c^VxB-- j 
ot eo 



'dt 



jp ) , 



= -V X E, 



(4.2) 
(4.3) 



where je and jp are the current densities of the electrons 
and protons. 

In our 2D simulations the y dimension is chosen to be 
invariant. The system is periodic in the z direction, so 
that a particle that leaves the system through the right 
or left boundary (see Fig. [1]) appears on the opposite 
boundary. Thus, it is not necessary to make the system 
very long in order to handle the whole particle trajec- 
tory from entering to ejection. The current sheet half- 
thickness, Lx, is equal to 1 m, while the width of the 



whole simulation region along x is chosen to be 20 m, 
in order to avoid any influence of the boundaries on the 
particles inside the RCS. Plasma is continuously injected 
from the x — ±10 m sides of the simulation region with 
the rate nEyo/ B^q. 

In order to avoid numerical instabilities in the PIC 
method, the following constrains need to be satisfied: 



cAt < A^, 
At < 0.2a;-\ 
AC < Xd, 



(4.4) 
(4.5) 
(4.6) 



where At is the time step, A^ is the grid step in any 
direction, c is the speed of light, — [Airne^ /rrig)^/'^ is 
the electron plasma frequency and Ad = (fcr/47rne^)^/^ 
is the Debye length. To satisfy these conditions while 
keeping the code running time at a reasonable level, we 
use a reduced plasma number density n = 10^° m"'^ in 
our simulations. Also the proton-to-electron mass ratio 
is reduced to nip/me = 100 in order to keep the proton 
acceleration time within reasonable computational limits. 
The spatial simulation grid has from 10 to 100 cells in the 
z direction and 100 cells in the x direction with Az = Ad, 
Ax = Ad/5 and 100 superparticles per cell in average, 
while each superparticle represents 2 ■ 10'' real particles. 
The time step is 6 • 10"^*^ s. 

As it was mentioned in Sec.[Tl we assume the existence 
of a background plasma with the real coronal density of 
10^^ m"'^, and which is not a part of PIC simulation, 
due to the limits of a current power even in large com- 
puter clusters. This assumption is required to provide the 
background current jy — — l//io • dBz/dx to sustain the 
equilibrium with the background magnetic configuration 
given by Eqs. (|2.2ll2.4p . that will produce electron veloci- 
ties Vy — l/{iione) ■ dBz/dx. For the real coronal plasma 
density the electron velocity which forms the equilibrium 
background current is 5 ■ 10^ m/s. This velocity is a few 
orders of magnitude lower than the velocities gained by 
the electrons at acceleration in the drifted electric field, 
thus it will not change significantly the electron distribu- 
tions and is neglected in the current approach. By doing 
so we, possibly, eliminate only some low frequency in- 
stabilities (like ion-sound wave), which are not a part of 
our present study. However, this approach allows us to 
investigate the acceleration of a smaller number of par- 
ticles (n = 10^0 m"^) accelerated by the drifted electric 
field by taking into account the plasma feedback to these 
particles motion. 

The simulations are supposed to run until a quasi- 
stationary state is reached. This running time corre- 
sponds to the acceleration time of the slowest particles 
(protons) and is usually less than 10"^ s. Also, since 
we are interested in the stationary state, and the back- 
ground magnetic field in assumed to be in equilibrium 
with the the high density plasma, the initial state of the 
simulated plasma particles is not important. In practice, 
all particles that are initially present in the simulation 
region would be accelerated, ejected and replaced by the 
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FIG. 4: Snapshots of the plasma particles on the x-Vz phase plane (PIC simulations). Current sheet parameters are 

the same as in Fig. [2j 



injected particles before the stationary state is reached. 
Thus, to speed up the simulation we do not initially have 
any particles in the simulation region. The injected par- 
ticles are assumed to have a thermal distribution with 
the typical coronal temperature of 10^ K. 



4.2. Comparison with the test particle simulations 

The PIC simulations proved that the RCS magnetic 
configuration may accelerate plasma particles even if the 
plasma feedback is taken into account. The results of 
simulation obtained for an extremely low number density 
of particles (10^ m^'^) are plotted in Figs.l4alandl4bl This 
low density simulation is only performed to verify that 
the PIC code can repeat the results obtained in the TP 
simulations (see Fig. [5]). Indeed, the formation of beams 
of accelerated protons (Fig.|4a|) and electrons (Fig, lib)) is 
clearly seen. Similar to the TP simulations, the electrons 
and protons are ejected into the opposite semispaces with 
respect to the midplane {x = 0). 



The protons are unmagnetised inside the current sheet, 
so the ejection velocity of the "bounced" and "transit" 
protons are almost the same. On the other hand, the 
electrons are much more magnetised. Hence, similarly to 
what was shown in the TP simulations, the "bounced" 
electrons are not able to reach the midplane and, thus, 
they gain less energy than the "transit" electrons. The 
electric field caused by the charge separation is quite 
weak in this case, about 17 V/m. 

In our further simulations the plasma density is ac- 
cepted to be higher (10^° m"-^), but it is still several or- 
ders of magnitude lower than the typical coronal den- 
sity (10^^ — 10^^ m~'^), which is unreachable due to the 
computational limitations. The higher density simula- 
tions (Figs. [3c] and I4d[) reveal more differences between 
the TP and PIC approaches. Proton dynamics in PIC 
(Fig. He]) has not changed much because of the increased 
density, they are still ejected mainly into one semispace 
{x < 0) with respect to the midplane. The acceleration 
rate of protons coincides with the theoretical value given 
by Eq. (|3.ip . and their trajectories are close to those ob- 
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FIG. 5: The fraction of protons that are ejected to the 

X > semispace as a function of ByQ (dashed Hne, 
B^o = 2- 10-5 T) and (sohd hne, Byo = 10''^ T). 
All other parameters of the current sheet are the same 
as in Fig. [4cl 



FIG. 6: Electric field Ex induced by particles in the 
PIC simulations for different values of Bxo and ByQ. All 
other parameters of the current sheet are the same as in 
Fig.H^ 
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tained from the TP simulations (see Fig. [2]). However, 
the proton orbit at the acceleration phase is wider (~ 3 m 
versus ^ 1 m in Fig. [5]). Also, near the midplanc the pro- 
tons form the structure with a higher density on the sides 
and in the centre (Fig. |4c]). A similar structure in the 
proton density inside the RCS was obtained by Zharkova 
and Agapitov [18]. 

In addition there is a small number of protons that 
are ejected to the x > semispace (Fig. where elec- 
trons are normally ejected in TP simulations. Plot in 
Fig. [5] shows the fraction of the protons that are ejected 
to the a; > semispace as a function of B^o and ByQ. 
The dependence on ByQ qualitatively coincides with the 
asymmetry rate dependence obtained by Zharkova and 
Gordovskyy In particular, for a very small ByQ in 
both PIC and TP simulations the protons are ejected 
symmetrically with respect to the midplane. However, 
in PIC simulation all the protons are ejected to the one 
semispace {x < 0) when ByQ> 1 ■ 10"^ T « B^q, while 
Zharkova and Gordovskyy [9| have shown that in the 
TP simulation accelerated particles are fully separated if 
ByQ > 1.5 • 10~^BzQ. Thus, the particle trajectories have 
smaller asymmetry in PIC simulations compared to the 
TP ones. 

The electron dynamics in PIC simulations, as it is seen 
from Fig. l4dl is different from what was predicted by the 
TP simulations. Instead of a formation of the beams with 
a narrow energy distribution, the accelerated electrons 
are found to gain a wide range energy spectrum. More- 
over, the electrons become ejected mainly into the same 
side as protons. Also the energies of ejected electrons in 
PIC simulations are smaller than in the TP simulations. 
In order to understand these new effects, let us study the 
electron and proton dynamics in more details by recon- 
structing the trajectories of particles in the RCS. 
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FIG. 7: Charge density p as a function of x deduced 
from the PIC simulation. Current sheet parameters are 
the same as in Fig. He) 



4.3. Induced electric field and particle trajectories 

The differences between the TP and PIC simulations 
could only occur because of the additional electric and 
magnetic fields induced locally by the accelerated parti- 
cles. The PIC simulations has shown that the induced 
magnetic field B is much smaller than the background B. 
On the other hand, the induced electric field E is essen- 
tial. Its absolute value is even larger than the background 
(drifted) field Ey induced by a magnetic reconnection (see 
Eq. jMl)). 

Since our 2D system is invariant in the y dimension 
there is no charge separation in this direction. Thus the 
induced Ey can be caused only by the time variation 
of the magnetic field and is found to be small. In this 
subsection we consider the electric field E^ which is per- 
pendicular to the current sheet. Fig. [5] is a plot of E^ 
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(a) Two protons, injected at x = ±10 m with zero initial velocities. 




(b) Electron, injected at z = — 10 m with initial velocity 
14 = -5 ■ 10*^ m/s. 
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(c) Electron, injected at x = 10 m with initial velocity 
\4 = -1 ■ 10^ m/s. 
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(d) Electron, injected at z = — 10 m with initial velocity 
= -5 ■ 10« m/s. 



FIG. 8: Trajectories of particles on the Vz-x phase plane during acceleration inside a current sheet. The trajectories 
were obtained in the TP simulations, with the electric field Ex{x) taken from the PIC simulation (Fig. [6] solid line). 

Current sheet parameters are the same as in Fig. He] 



as a function of x averaged over the z coordinate for dif- 
ferent values of B^o and ByQ. This field appears due to 
separation of the electrons and protons across the cur- 
rent sheet, which leads to a local non-neutrality of the 
plasma. The field becomes stronger when B^o decreases 
or ByQ increases. Also, when Byo the electric field 
Ex and the system itself become symmetric. The dis- 
tribution of a charge density p{x), which generates the 
polarisation field E^ix), is shown in Fig. [T] Note, that 
the average charge density over the simulation region is 
close to zero, which mean that the current sheet as a 
whole remains electrically neutral. 

In order to reconstruct the particle trajectories, we use 
the TP code, where the induced electric field Ex obtained 
from PIC is added to the background electro-magnetic 
configuration given by Eqs. (|2.2j|2.5p . The trajectories 
of the two protons in the x-Vz phase plane is shown in 
Fig. [Sal These trajectories are very similar to those ob- 



tained in the TP simulations without E^ (see Fig. [2]). The 
only difference is that during the acceleration phase the 
"bounced" proton has a wider orbit. It is clear now that 
this wide orbit is responsible for the two smaller peaks 
at about ±1 m in the charge density plot (Fig. [7]). On 
the other hand, the narrow orbit of the "transit" proton 
forms the central peak in this plot. 

The trajectories of electrons are much more compli- 
cated. Firstly, let us consider an electron that enters 
from the x < semispace (Fig. l8b[) . The dynamics of this 
electron is similar to the dynamics of the "transit" elec- 
tron in Fig. [5]- the electron drifts towards the midplane, 
becomes accelerated and ejected to the x > semispace. 
However, the polarisation field Ex{x), which extends be- 
yond the current sheet and has a component parallel to 
the magnetic field, decelerates the ejected electron. For 
the chosen magnitudes of B^ and By the majority of 
electrons can not escape to the a; > semispace, instead 
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(c) B^o = 10-4 T, Byo = 10-4 T (d) 5,0 = 5 • lO-^^ T, Byo = 10-4 T 

FIG. 9: The energy spectra of the ejected electrons. Dashed hne ~ Eq. (|3.ip . dotted hne - Eq. (|3.2p . dash-dotted hne 
- mean energy of the ejected electrons. Current sheet parameters are the same as in Fig. |4cl 



they are dragged back towards the current sheet and be- 
come indistinguishable from the electrons entered from 
the a; > semispace. 



The electrons that come from the positive x side 
demonstrate a rather different dynamics in comparison 
with the case when = 0. It turns out that the electron, 
which is "bounced" from RCS in the absence of Ex, can 
now reach its midplane. In the vicinity of the midplane, 
the electron becomes unmagnetised and oscillates with 
the gyrofrequency determined by By (Fig. [8c]) and, after 
some time of oscillation, the electron is ejected. If the 
electron initial velocity is small, it can be quasi-trapped 
inside the RCS (Fig. I8dp . Such electron is accelerated on 
the midplane, ejected from it, then decelerated outside 
the RCS and returns back to the midplane. This cycle 
is repeated for several times, until it finally gains enough 
energy to escape the RCS. Since the magnitude of the 
polarisation field E^ (x) is smaller at x < than at a; > 
(see Fig. [5]), it is easier for the electron to escape to the 
a; < semispace. Thus, most of the electrons are ejected 
to the same semispace as protons. 



4.4. Energy spectra 

The energy distributions of accelerated electrons for 
different magnitudes of transverse, B^o, and guiding, 
ByQ, magnetic components are calculated for those 
electrons outside the RCS in the negative semispace 
(—10 m < a; < — 5 m) and plotted in Fig. [9l Only escap- 
ing electrons are taken into account, i.e. the electrons for 
which the velocity component along the magnetic field 
projected on the x axis is larger than the E x B drift 
velocity projected on the x axis. For the accepted con- 
figuration and in the limit of a small B^ , this corresponds 
to the electrons with Vz > Eyo/Bxo- 

It can be seen from Fig.[n]that for any B^o and Byo the 
ejected electrons form a wide single-peak energy distribu- 
tion with a width of the order of the mean energy. This is 
contrary to the TP simulations and the lower density PIC 
simulations, where the two narrow energy electron beams 
are formed (see Fig. [2] for the TP simulation and Fig. I4bl 
for the lower density PIC simulation). The presence of 
the polarisation electric field shifts the high energy peak 
towards the low energy one, while the latter is expanded. 

If the guiding field. By, is negligibly small (Fig. [9a|) . 
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FIG. 10: Electric field Ez induced by particles in PIC 
simulation {B^q = 10"^ T, Byo = 10"* T, 
S^o = 4 • IQ-^ T, EyQ = 250 V/m, mp/m^ = 10, 
n = 10^2 ^-3y 
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FIG. 11: Distribution function of electrons calculated 
for the region indicated by the dashed rectangle in 

Fig.iini 



then the mean energy of the ejected electrons coincides 
with the analytical value given by Eq. p.ip . For the 
stronger guiding field (Fig. I9bp the mean energy is some- 
what higher than the energy of the slow beam in the TP 
simulation given by Eq. (|3.ip . and substantially lower 
than the energy of the fast beam given by Eq. p.2p . 

The energy spectra dependence on is shown in 
Figs. [9c] and I9dl They indicate that the mean energy 
strongly depends on B^ as it is predicted by Eqs. p.ip 
and (j3.2p . The lower B^ the higher their mean energy 
and the wider the energy distribution. 



4.5. Generation of the Langmuir waves 

Another interesting effect observed in the PIC simu- 
lations is the excitation of waves induced by accelerated 
particles. As it can be seen from Fig. [TOl the com- 
ponent of the induced electric field is structured with a 
characteristic length scale of about Xwave ~ 2 m. This 
structure propagates in time in the positive direction 
of the z axis. The speed of this propagation is about 
Vwave ~ 1-3 ■ 10^ m/s, which makes the temporal pe- 
riod of oscillations about T^ave ~ 1.5 • lO"*" s. Since the 
plasma frequency for n — 10^^ m^"^ is 9.1 • 10^ Hz, the 
generated wave is the Langmuir wave. Note also that 
the oscillating component of the excited wave E^ is par- 
allel to the direction of propagation, which corresponds 
to the polarisation of the Langmuir wave. 

In order to clarify the mechanism of instability, let us 
plot a Vz distribution function of electrons in the region 
where the instability occurs (this region is indicated by 
the dashed rectangle in Fig. fTU)) . Fig. [TT] indicates that 
the electrons have the typical "bump-in-tail" unstable 
distribution. The range of velocities Vz for which the 
derivative df jdVz is positive is from 1.3-10^ to 2-10^ ni/s, 
these values correspond to the phase velocity of the Lang- 
muir wave Vwave deduced earlier. 

The generation of the Langmuir waves is not the only 
possible instability in an RCS. For example, Drake et al. 
[23| have shown the possibility of whistles waves excita- 
tion in thin current layers. However, in our simulations 
we can observe only waves with length scales less than 
the size of the simulation domain. Also, since the current 
simulation system is 2D in space (and 3D in velocities), 
the waves which are essentially 3D cannot be generated. 
The investigation of other instabilities, apart from Lang- 
muir waves, will be a subject of the forthcoming study. 

5. CONCLUSIONS AND DISCUSSIONS 

In this paper we present the PIC simulations combined 
with the TP approach for particle acceleration in a re- 
connecting current sheet with given magnetic field topol- 
ogy. The acceleration is studied in a small part of the 
current sheet with a 3D background magnetic field and 
a drifted electric field induced by the magnetic diffusion. 
The simulations are carried out for thin RCS of 2 m. The 
main (reconnecting) component of the magnetic field is 
10^'^ T, while the other two components are varied in 
wide range. 

The background magnetic and electric fields are con- 
sidered to be stationary while varying across the current 
sheet and being constant in other dimensions. The PIC 
simulations is used to study particle acceleration in the 
combined background electro-magnetic fields and those 
fields induced by the accelerated particles. The TP ap- 
proach is used to study how these additional fields affect 
the particle trajectories inside the 3D RCS, that, in turn, 
helps to understand the PIC results. 
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The PIC simulations carried out for a lower plasma 
density (10'' m~^) reveal the same particle trajectories 
as in the TP approach. In particular, for the RCS with 
lower density we confirm the separation, or asymmetry, 
of particle beams with the opposite charges, which are 
ejected into the opposite semispaces with respect to an 
RCS midplane. For a higher plasma density (10^'^ ni~'^), 
the accelerated particles are found to produce the addi- 
tional electric fields: a polarisation field, E^, caused by 
the charge separation of accelerated particles across the 
current sheet and by a turbulent field of Langmuir waves, 
Ez, caused by the "bump-in-tail" energy spectra of the 
accelerated electrons. 

The magnitude of the polarisation field Ex is shown 
to be very high that leads to essential modification of 
the trajectories of electrons, and, to some extent, those 
of protons. There is a difference in the acceleration of 
the "transit" particles (those which enter from and are 
ejected to the opposite semispaces towards the midplane) 
and " bounced" particles (those which enter from and are 
ejected to the same semispace towards the midplane). 

In particular, during their acceleration, the "bounced" 
protons are found to have wider orbits than the "transit" 
ones. It was also found in PIC simulations that the pro- 
tons are ejected less asymmetrically with respect to the 
midplane compared to the trajectories found in the TP 
simulations. For example, in the TP approach the pro- 
ton ejection is fully asymmetric (all protons are ejected 
to the same semispace) if the guiding field ByQ is larger 
than 1.5 • IO^^Bzq, while in the PIC simulations this crit- 
ical guiding field is of the same order of magnitude of the 
main component Bzq. 

In the current RCS topology, if there is an asymmetry 
of the accelerated protons at ejection, the accelerated 
electrons in the PIC simulations are shown to be ejected 
to the same semispace as the protons, i.e. there is no 
non-local charge separation as predicted by the TP sim- 
ulations jSi]. The electrons, which are "bounced" in the 
TP approach, are in the PIC simulations dragged by the 
polarisation field E^ to the midplane. At the same time, 
the "transit" electrons, which are ejected to the opposite 
semispace from protons in the TP approach, cannot leave 
the simulation region in PIC, i.e. they become dragged 
by the polarisation field E^ back to the RCS. Since in PIC 
the electrons are eventually ejected into the same semis- 
pace as protons, the charge separation can only happen 
locally, about 10 m around the midplane. This is valid 
for the small simulation region considered in the current 



paper, while the PIC simulations for the larger region 
could result in separation of high energy particles across 
the RCS. 

The energy distribution of ejected electrons in the PIC 
simulations is found to be rather wide, i.e. the width 
of the distribution is of the order of magnitude of the 
mean energy, which is different from the energy spectra 
found in TP approach [9]. However, the mean energy of 
electrons is found to be consistent with the TP simulation 
and the analytical estimations for a weak guiding field. In 
spite of the strong guiding field can significantly enhance 
the ejection energy of the "transit" electrons in the TP 
simulations, in the PIC simulations this enhancement is 
shown to be much less pronounced if the polarisation field 
Ex is taken into account. 

PIC simulations reveal the additional periodic electric 
field Ez- This field is found to be formed by the waves 
generated by the beam of accelerated electrons with 
the unstable "bump-in-tail" velocity distribution. For 
the magnitudes of the background electric and magnetic 
fields considered in the current study, the phase speed 
of generated waves is w 1.3 • 10^ rn/s, the wave length 
is w 2 m and the period of oscillations is « 1.5 ■ 10^^ s, 
which corresponds to the Langmuir wave frequency. The 
turbulent electric field Ez of this waves is likely to be 
the reason of the wide energy distribution of the ejected 
electrons derived from the PIC simulation. 

Therefore, we show that the feedback of the ambient 
plasma in an RCS in a form of electric fields produced 
by the accelerated particles is rather strong and modifies 
substantially the particle trajectories and energy spectra 
gained during their acceleration. These effects are essen- 
tial even for a small fraction of the particles simulated 
by PIC in the present study. The full particle simula- 
tion for a real plasma density can reveal further effects 
in the plasma dynamics during a magnetic reconnection, 
that will be a scope of the forthcoming studies with the 
increased computer power. 
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